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One gluon, two gluon: multigluon production via high energy 
evolution 
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Abstract: We develop an approach for calculating the inclusive multigluon production within the 
JIMWLK high energy evolution. We give a formal expression of multigluon cross section in terms of 
a generating functional for arbitrary number of gluons n. In the dipole limit the expression simplifies 
dramatically. We recover the previously known results for single and double gluon inclusive cross 
section and generalize those for arbitrary multigluon amplitude in terms of Feynman diagramms of 
Pomeron - like objects coupled to external rapidity dependent field s{rj). We confirm the conclusion 
that the AGK cutting rules in general are violated in multigluon production. However we present an 
argument to the effect that for doubly inclusive cross section the AGK diagramms give the leading 
contribution at high energy, while genuine violation only occurs for triple and higher inclusive gluon 
production. We discuss some general properties of our expressions and suggest a line of argument to 
simplify the approach further. 
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1. Introduction 



The study of high energy hadronic scattering and the physics of saturation has been a very active 
area of research in recent years [|l|, ^, ^, |5|, ^ |^ . The effort has bifurcated recently with one branch 
concentrating on understanding dynamics of Pomeron loops |9|, [l^, |ll|, |l2|, |l^ while the other 
one pushing towards better understanding of more detailed aspects of JIMWLK evolutionj^, ^, such 
as relation to the QCD Regge theory ^ and calculation of observables less inclusive than the 



total cross section [|16|, |17|, |18|, |19[. The latter line of research seems especially relevant since we 
understand by now that the total cross section cannot be calculated perturbatively at asymptotically 
high energies, even if the Pomeron loops are taken into account p^ . 

The present paper is devoted to a detailed discussion of multigluon production in the JIMWLK 
approach. The study of gluon production at high energy was initiated in a paper by Kovchegov 
and Tuchin |21| which calculated the single gluon inclusive cross section within the dipole model 



approximations. This was later extended beyond the dipole model [p2[,[17|. The expression for the 



double gluon inclusive production in the dipole limit with an arbitrary rapidity interval between the 
two gluons was obtained by a real tour de force calculation in [l^ , while a simpler calculation of two 



gluon production at the same rapidity but beyond the large Nc limit was performed in |17]. The 
double gluon inclusive cross section was also calculated by Braun |^] in the framework of the BFKL 
Pomeron [23| calculus. One should note that the results of and |24| are different, as the AGK 



cutting rules |25] are violated in |16| while |24] suggests that the AGK violating contributions are 
small at high energy. 

All the above calculations have been performed on a case by case basis, and at the moment 
there is no unified approach that would streamline calculations of an arbitrary observable of this 
type. Recently [l^ we have introduced such a unified formalism for calculation of semi-inclusive 
observables such as various diffractive and elastic cross sections. In the present paper our goal is to 
extend [19| to calculation of multigluon production. 



The structure of this paper is the following. In Section 2 we recall the general formalism of the 
high energy evolution as it arises from the evolution of hadronic wave function and applies to semi- 
inclusive observables. In Section 3 we extend the formalism to include multigluon production. We 
write down a formal expression for the generating functional which generates the inclusive probability 
of production of n gluons at arbitrary rapidities Yi and arbitrary transverse momenta ki. We discuss 
some properties of this expression and show that it reproduces the known result for the single gluon 
inclusive case. 

In Section 4 we specify to the dipole model limit in which the problem simplifies considerably. 
In this limit we can develop a simple diagrammatic technique in terms of the propagators of BFKL 
Pomeron-like objects propagating in the external field of the dipole s-matrix s{r]). We confirm that 
for double inclusive gluon production our result agrees with |16|. We stress that not only the BFKL 



Pomeron but also other similar objects propagate in the relevant diagrams if the number of produced 
gluons is greater than one. It is the presence of these additional objects that violates the AGK 
cutting rules. The actual calculation of the multigluon cross section requires the knowledge of the 
propagators of these "cousins" of the BFKL Pomeron. We find the asymptotic behavior of these 
propagators for small and large separations and show that the crossover must occur on the scale of 
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the inverse saturation momentum of the target. Based on the behavior of the propagator and of the 
gluon emission vertices we argue that the leading high energy contribution to double inclusive gluon 
production still preserves the AGK cutting rules thus reconciling the results of ||l^ and [24|. We 
show however that for triple and higher gluon production the AGK rules are genuinely violated. We 
formulate an ansatz which incorporates the known asymptotics of the relevant propagators with the 
correct crossover scale which we believe should be a good approximation to physics. We hope that 
numerical calculations using this ansatz may be feasible. 



2. High energy scattering: the general setup 
2.1 Evolution of a wavefunction 

In this section we recall the high energy evolution of the hadronic wave function as discussed in 

iO, MM- 

Consider an arbitrary projectile whose wave function in the gluon Fock space can be written as 

|*)yo = ^[af(:r)]|0)y„. (2.1) 

The gluon creation operator depends on the transverse coordinate, and also on the longitudinal 
momentum fc"*". The projectile is assumed to have large energy. The gluon operators which enter 
eq.( |2.1D therefore all have longitudinal momenta above some cutoff A. We will refer to these degrees 
of freedom as "valence", |^') = \v). Henceforth we omit the dependence on longitudinal momentum 
in our expressions, as the momentum enters only as a spectator variable and only determines the 
total phase space available for the evolution. The wave function in eq. (|2.l| ) is arbitrary. 

When boosted, the valence wave function gets dressed by an evolution operator 17 ^] , which 
depends on the valence degrees of freedom {p) as well as on newly created soft modes (a). The 
evolved wave function has the following structure 

= nyip, a) \v)yo ; \v)yo = \v) ^ K) (2.2) 

where |0a) is the Fock vacuum of the Hilbert space of the soft modes. The evolution operator Q is 
explicitly known in the dilute limit {p — 0) only. In this limit, = C, where C is the coherent 
operator 



Cy = Expji I d2^6^«(x)^'^^^-^^^^ [<(A:+,x) + (fc+,x))] | . (2.3) 

The creation operators a'^{k^) create gluons with soft momenta, which are not present in the valence 
state \v). The field b depends only on the valence degrees of freedom. It is determined as the solution 
of the "classical" equation of motion 



dib^ + ge''''^b°{x)bf{x) = gp''{x) 
b] - djb^ + ge''H\b] 



[dM - djb^ + ge''H\b'^] = (2.4) 
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where p'^{x) is the color charge density operator of the projectile. The coherent operator C dresses 
the valence wave function by a cloud of the Weiszacker- Williams gluons. Operatorially 

C^A'i{k+,x)C = AUk+,x) + ^bUx). (2.5) 

The evolution of a generic non-diagonal matrix element of any observable O is given by 

d{v\d\v') .. {v\^}\.d{p + Sp)nY\v') - {v\d(.p)\v') i\m^n( U'\ cya^ 

Here both states \v) and \v') are boosted by the same amount. With the explicit form of the evolution 
operator C the evolution Hamiltonian is obtained as 

A] = j^^Qt{z)Qnz) (2.7) 

where the one gluon emission amplitudes Qf{z) is defined as 

QUz) = R^\z)bUz) - bUz), (2.8) 

with the "dual Wilson line" R 



Riz 



ab 

,ab 



Pexp / dz'T" 
Jo 



lo ^P'^iz, z ) 



(2.9) 



Here bL,R are solutions of the classical equations of motion with sources given by the generators of 
the left/right color rotations Jl,r defined as 

4(x) = -tr |i?(x)T"^^| , Jlix) = -tr |™(x)^^| , Jl{x) = [R{x) Jr(x)]« . 

(2.10) 

The rotations have an alternative representation in terms of the charge density operators: 

4(x) = p-ix,x~ = iy, J£(x) = p'^{x,x'=0). (2.11) 

The KLWMIJ+ Hamiltonian reduces to the KLWMIJ Hamiltonian when we keep the linear approx- 
imation for the classical field b: 

6r(.) = |-|d\|^p». (2.12) 

and then 

Qt{z) = J d'x^^lR'^^z) - R'^\x)]jUx). (2.13) 

It is the last form that we will use explicitly throughout this paper although generalization to KL- 
WMIJ+ does not present any difficulties. 

Within the more general approach including the Pomeron loops, the evolution Hamiltonian can 
be written in terms of the n-gluon "emission amplitudes" Qn []T5|. 
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C 7 

Hip^y] = — m^\.{p,a)^Y{Rp,Ra)MY=Y, = Qi[R]Qn[R] ■ (2.14) 

Each Qn depends on transverse coordinates and has a set of n Lorentz and color indices which we 
do not indicate exphcitly. Unfortunately the explicit form of Q„ is not known at present. 

Throughout this paper we will be thinking of the diagonal matrix element as being 

computed through the path integral representation 

{v\d\v) = [ Dpd W[p]. (2.15) 



As discussed in ||9|, 13] this representation requires endowing p(x) by additional ordering coordinate, 
but this subtlety is of no concern to us in this paper. 

2.2 High energy scattering 

Let us now consider the process whereby the projectile discussed above scatters on a target. The 
target is described by some distribution of target chromoelectric fields. We are working in the 
lightcone gauge natural to the projectile wave function, A~ = 0. In this gauge the high energy 
scattering matrix of a single gluon projectile at transverse position x on the target is given by the 
eikonal factor 



S{x;0,x-) = V exp{i / dy- T" aUx,y-)} ; S{x) = 5(x;0,l) 

Jo 



(2.16) 



The field ot is the large A'^ component created by the target color charge density. It obeys 
the classical equation of motion and is determined by the color charge density of the target pt{x) 
vial2|, I, 

a^(x,x-)r'^ =g^ ^{x-y) {5t(y; 0, x") p^(y, x") T'^ 5(y; 0, x")} . (2.17) 

For a composite projectile which has some distribution of gluons in its wave function the eikonal 
S-matrix can be written in the form analogous to ^(x), see |10] 



Y.PP[aT] = {P\S\P) = j dppW^'lpp] exp!^ij\y- j d^x p''p{x,y-) aUx,y-)^ (2.18) 

where S is the second quantized 5-matrix operator of QCD. The quantity pp{xi) is the color charge 
density in the projectile wave function at a given transverse position, while l^''^[p] can be thought 
of as the weight functional which determines the probability density to find a given configuration of 
color charges in the projectile. For a single gluon p"(xj) = r"(5^(xj — x^), and eg. ( p. 18 ) reduces to 
eq. (ra . 

The total 5-matrix of the scattering process at a given rapidity Y is given by 

S{Y) = y^a^ VFfjaT(x,x-)] S^^yjar(x,x-)]. (2.19) 
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In Eq. ( 2.1S| ) we have chosen the frame where the target has rapidity Yq while the projectile carries 
the rest of the total rapidity Y . Lorentz invariance requires S to be independent of Yq. 

The high energy evolution of the projectile wave function discussed in the previous subsection 
leads to the following expression for the evolution of the forward scattering matrix element: 



dY 



j Da^ W^^[aTix,x-)] H 



S 



5 ax 



^pp 



(2.20) 



Here the (Hermitian) high energy evolution Hamiltonian H can be viewed as acting either to the 
right or to the left: 



d_ 
dY 



H 



ax, 



6aT 



^pp 



[«t] ; ^ W 



H 



6aT 



(2.21) 



As was shown in [IC] in order for the total S-matrix to be Lorentz invariant and symmetric between 
the projectile and the target, the evolution Hamiltonian H must be self dual. This is to say that the 
full evolution Hamiltonian which includes Pomeron loops has to be invariant under the Dense-Dilute 
Duality transformation 



a"-{x,x ) 



Sp"'{x,x )' 



(5q"(x,x ) 



-ip''{x,x ), 



S 



(2.22) 



However if one considers the situation where the target is large and the projectile is small (the dense- 
dilute limit), the symmetry between the target and the projectile is irrelevant. In the limit when the 
color charge density of the target is parametrically large {p'^ = 0(1 /ag)) the Hamiltonian is given by 
the JIMWLK-F expression |, 0] 



H 



JIMWLK^ 



Qt{^)Qi{z) 



(2.23) 



with Q[S] obtained from Q[R\ (|2.8D by the substitution R ^ S. In particular in the JIMWLK limit 



d^x [S''\z) - S'^'^x)] 4[S, x] . 



{x - zf 



(2.24) 



Thus for processes involving small perturbative projectile, the evolution of its wave function VF[pp] 
is governed by jj^^^^^^J while the evolution of the projectile averaged 5-matrix T,[S] is given by 

jjJIMWLK 



2.3 Semi-inclusive reactions 

Consider an observable O which depends only on the projectile degrees of freedom, which we want 
to calculate in the asymptotic state long time after the scattering has taken place. We assume that 
the target has been boosted to rapidity Yq and that the operator O depends on degrees of freedom 



between the rapidity Yq -\- y and Y. As discussed in |19] such an observable is given by 



id) = {TY,{PY-Yo-y\Cl{l - S^)Cy a Ct (1 - S)Cy\Py^y,^y)TY,) (2.25) 
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where \Tyq) is the wave function of the target boosted to rapidity Yq and \PY-Yo-y) is the wave 
function of the projectile boosted to rapidity Y — YQ — y. To calculate the evolution of this observable 
with respect to y at fixed Yq and Y — Yq — we consider first the same expression before taking its 
average over the target wave function 



0[S,S] = {PY^Y,-y\Cl{l - S^)Cy O Ct (1 - S)Cy\PY-Yo-y) 

The two are related by 



(O) = {TyMS,S]\s = s\Ty,) 



DSDSWl[S]5{S - S) 0[S,S]. 



The derivation of the evolution is similar to that given in the first subsection: 



dO[S, S] 



Mm 

dy A?/^o Ay 

We write Cy+^y = C^y Cy and expand Caj/ up to quadratic order in a: 



Oy+/\y[S,S] — Oy[S,S] 



CAy 



I + i I d^xh1{x) 



dk^ 



J d'xb'}ix)l^ 



A 7rl/2|A;+|l/2 

dk+ 



at{k+,x) + af(A;+,x)) 



leVA ^l/2|A;+|l/2 

We now recall the identities derived in 1 11 111: 



^{k+,x) + al^{k+,x)) 



(VI /St = Jt[S]{^P\S^ 

Combining the expansion (I2.29D and the identities (^ ) we obtain 

d 



Sp^l^) = Jt[S]S\^P); 
(VI = J^[5](V|5t. 



dy 



H3[S,S] = 
Let us introduce the notations 

Hi[S] = H-^IMWLK^g^ 



0[S,S] =-H3[S,S] 0[S,S] 



H2[S,S] = Hi[S] + Hi[S]. 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 

(2.31) 
(2.32) 



The Hamiltonian was first introduced in ||18[| . It is straightforward to prove that the Hamiltonian 
H3 obeys two useful identities. First, when acting on a function F which depends on the product 
S St only it reduces to the action of the Hamiltonian Hi : 



H3[S,S] F[SS^ = Hi[SS^ F[SS^]. 



(2.33) 



Second, when H3 acts on an arbitrary function ^[5*, S*] and then S is set to be equal S it is again 
equivalent to the action of Hi : 



H3[S,S] F[S,S^] 



s=s 



Hi[S] F[S,S^ 



(2.34) 
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Eq. ( 2.34 ) can be written in the following form 

J DSdiS - S)H3[S,S] F[S,S^ = Hi[S] j DS6{S - S)F[S,S^] (2.35) 



Referring back to eq.( 2.31| ) we see that for any observable the evolution operator between rapidity 
Yi and Y2 is 

UsiYi, Y2) = exp[-Hs{Y2 - Yi)] . (2.36) 

The preceeding discussion is easily generalized by considering a product of n + 1 local in rapidity 
operators Oi each one depending only on degrees of freedom at fixed rapidity Yq + Yi. We will call 
Oq (which plays a special role) the operator which depends only on the valence degrees of freedom 
of the projectile at rapidity Y. The expectation value for such an observable reads 

{On... do) = J DSDSWl[S]d{S - S) X 

X UsiYo, Yo + y„) On[S, S] UsiYn, y„-i) On-i[S, S]... Us{Yo + Yi,Y) Oo[S, S] . 

(2.37) 

The operator Oq determines what kind of observation is made on the valence degrees of freedom 
of the projectile. In the present paper we are interested in observables which are inclusive over the 
valence projectile degrees of freedom. For these observables (see ]l9[) 



Oo [s,s] = 1 - ^ris^ - ^ns] + ^ns^s] (2.38) 

and this is the expression we use in the rest of this paper. If on the other hand the process of interest 
is elastic over the projectile, we would have 

Ooeiastic[S, S] = {1- (1 - ^nS]) ■ (2.39) 

Finally for example jets produced in association with diffractive vector meson production at projec- 
tile's rapidity one has = |^') (^'| with ^' be a wavefunction of a final vector meson state. 



3. Inclusive gluon production 

3.1 Single inclusive gluon production 

Following the derivation of the previous section we define the operator Og corresponding to an 
inclusive production of a gluon with momentum ki at rapidity Ig + Yi 

Og[S,S] = A(o„|C,(l - 5t)Ct niki,y) Cy {1 - Cl\Oa)\y=Yo+n (3-1) 
u y 

h{k^,y) = dk+ a\\ki,k+)a1{ki,k+). (3.2) 

Je(n)+i'i) A 
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In the dilute limit with the operator C given in eq.(2^) the averaging over the state |0a) can be 
performed explicitly. This requires expanding the exponential operator C to second order in Taylor 
expansion. The result is 



(Pz (Pz 



Jki (z-z) 



Qt{z,[S])Qt{z,[S]) 



(3.3) 



with the gluon production amplitude Q the same as in eq.( 2.24 ). Substituting Eq. (3^) into equation 
Eq. ( 2.37| ) and remembering to set S = S in the end of calculation we obtain for single inclusive 
gluon production 



da 



dYi dkl 



(3.4) 



j DSDSW^^[S] 5{S - S) U^iYo, Yq + Y^) Og{ki) U^iYo + Yi,Y) £^[5^5] 



as only the last term in eq.(2.38) gives a nonvanishing contribution. Using the identity Eq. ( 2.34| ) we 
obtain 



da 



dYi dkf 



DSWl[S]Ui{Yo, Yo + Yi) J DS6iS - S) Ogiki)U3iYo + Y,,Y)i:^''[S^ S] 
DSDSWI^yAS] 6{S - S) Og{ki) U3{Yo + Yi,Y) E^i^^^] . (3.5) 



Using the fact that T,^^ evolves according to JIMWLK evolution we can write it as 
da 



dYi dkl 



I DSDSW^^+YAS]6iS - S)Og{k,)^^^y^_Y,[S^S]. 



(3.6) 



Substituting eq.(^^) and assuming color neutrality of the target this becomes (Iq = + 



da 



dYidkf 4 7r3 



ikiz-z) f {z-y)i 

{z — xy {z — nY 



DSWUS] tr \{Sl - Sl){S, - Sy) 



s=s 



(3.7) 



In order to make contact with the known results |21, 17] and also to streamline the future discussion 
let us change variables from S io U defined as: 



S = SU 



U = S: 



U^U = I. 



(3.8) 



In the end of the computation we will have to set U = 1. 
It is straightforward to derive the following identities 



Jr[S]\s,s = Jm\s,u; 

Jr[S]\s,s = Jr[S]\s,u + Jm\s,u; 

Using those we obtain 
n^(x, y;Y- Yo) ^ Ij^S, x] J^S, y] ^^"ly^ [S^ S] 



Jl[S]\s,s = S^' Jl[U]\s,u ; 

Jl[S]\s,s = Jl[S]\s,u + S'^' Ji[U]\s,u ■ (3.9) 



s=s 



jR[U,x]JUU,y]K-YAU] ^ . (3.10) 



U=l 
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The function obeys the BFKL equation. To see this one can use the identities ( 2.3C ) and express 
as a correlator of projectile color charge density operators 



n^{x, y) 



DppW^ipp] p%{x)p%{y) 



(3.11) 



In the dilute limit evolves with the KLWMIJ Hamiltonian, which when integrated by parts to 
act on pp reduces to the BFKL equation. The final result for single gluon inclusive production reads 

_ as f ^ik(z-z) f iz-x)i {z-y)i 
dyi dkl vr J^ ^ J^y {z - x)^ {z - yf 

X nP{x,y-Y-%) [{T.^yh, + {T^M 
with (T) denoting an ^-matrix of a gluonic dipole: 

{T.,y)Y, ^ j DSW'y^[S]tT[SlSy\ 

3.2 Multi-gluon production 



(3.12) 



(3.13) 



The formal generalization of single gluon inclusive cross section to an arbitrary number of gluons is 
straightforward. We first introduce a generating functional Z: 



I DSWliS] I DS6{S - S) Ve- '^"[^^"^^ o,{k)j{n,k)] ^PP^g\g^ (3 



The n-gluon cross section is calculated as ^ 

da 



j=0 



(3.15) 



dYi dkl... dYn dkl ^j{Yl , ^1) ^3{yn, kn) 

It is interesting to note that the generating functional Z can be also used to define diffractive 
processes with rapidity gap(s). To impose a rapidity within an interval between Yi and I2 we 
simply set ^{7), k) = 2 within this interval and zero elsewhere. This procedure changes the evolution 
Hamiltonian within this gap to i?2- As discussed in [^] this is precisely what is needed to forbid 
gluon emission within the rapidity gap. 

The formal expression eq.( 3.15| ) is unfortunately of little help as it stands. We can make a 
little bit of headway in understanding its structure by introducing the variable U as in the previous 
subsection. A straightforward calculation allows one to express both the Hamiltonian and the gluon 
production operator as follows 



H3 = H,[S] + H,[U] + Vsu[S,U] ; Vsu[S,U] = 2 [ Qt[S]Qt[U] , 

J z 



(3.16) 



Oy{k) = -f 



x,y,z,z 



'"^'-'^J^x,z,y,z\jt[x,U] \{Sl - Sl){S,U: 



Sy Uy] 



ab 



+ J^[x, S] {Si - Si) (S, U, - Sy Uy) J'j,[y, U] 



(3.17) 



^Note that instead of introducing the current j we could introduce two different currents coupled separately to Q[S] 
and Q[S]. That would be in analogy with the Schwinger-Keldish formalism discussed recently in |29[. 
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where the kernel N is defined in eq.( 4.10| ). 

In the new variables the generating functional Z reads: 



DSDUWy^\S]5{U - l)Ve' 



(3.18) 



and the multigluon production cross section 
da 



dYi dkf ... dYn dkl 



- = I DSDUWl^^y^[S]5{U - 1) Og{ki,[S,U])U^i{Y2, li) x 

n J 



X Og{k2AS,U])U^{Y^, Y2)...Og{knAS,U])T.^1y^^y^p] 



(3.19) 



Here we have used the fact that both, the constant function and the delta function are eigenfunctions 
of the operators Q"[x, S\ with vanishing eigenvalue. That allowed us to drop the [/-dependent terms 
in leftmost evolution operator and the 5-dependent terms in the rightmost evolution operator, and 
thus simply to evolve VF"^[S'] and T:^^\U] to the rapidity of the closest emitted gluon. A further 
slight simplification that follows from the same argument is that in Og{kn) we can drop the second 
term which involves /LI'S], while in Og{ki) we can set U = 1; 



vr 



x,y,z,z 



*'"(^-')A4,.,,,.- Jl[x,U] {Si - Sl){S,U, - SyUy) Ji[y,U] (3.20) 



ab 



o 



Jx. 



'^(^-^")A4,.,,,.- Jl[x,U] {Si - Sl){S, - Sy 4[y,C/] + 



ab 



+ rn[x,S] {Si - Sl){S, - Sy 



ab 



(3.21) 



We note the structure of eq.( |3.19| ) is rather curious from the point of view of the Hilbert space of the 
JIMWLK Hamiltonian [|l5| . The expression involves the off diagonal matrix element of the evolution 
operator with insertions of Og's between two states which are products of factors depending on S and 
U: Wf^+yJ5] 6{U - 1) and ^^^Y„_y,P]. The factor 6{U - 1) on the left is nothing but the Yang 
ground state of H-^^^^^^^[U], while the unity on the right is the Yin ground state of H-^^^^^^^[S] 
1 15]. Thus the expression for the multigluon production incorporates some elements of the Dense- 
Dilute Duality symmetry. This can possibly give a clue for future progress in calculating multigluon 
production. 

At the moment however despite the simplification of eqs.( p.2d ,3.21) the calculation of multigluon 
production within the general JIMWLK looks like a very difficult problem. The Hamiltonian H-^ 
couples S and U in a nonperturbative way and this coupled problem has to be solved. We note that 
since the evolution of U is eventually projected onto the Yang state, one can develop perturbation 
theory for the U part of the evolution. In particular writing 



U{x) = exp{i (3" T"} 



(3.22) 



and considering /? to be small we can perform expansion in powers of /?. This is useful since the 
Yang projection sets /3 to be equal to zero at the end of the calculation. As we shall see in the next 
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section this strategy is very useful within the dipole model approximation. Here however, even if /? 
is considered to be small, its interaction with S in H3 is not suppressed. Thus the first step in this 
perturbation theory requires solution of the coupled S and /? theory. This problem is beyond our 
abilities at the moment. We can make however quite a bit of headway by restricting to the dipole 
limit. This we do in the next section in the hope to come back to the general problem in the future. 



4. Multi-gluon production in the dipole limit 



In this section we discuss the dipole limit of the formulae derived in the previous section. We will see 
that in this limit the multigluon production amplitudes simplify dramatically. The starting point of 



our discussion is eq. ( 3.14 ). 
4.1 The dipole limit 

The dipole operator for S and S are defined conventionally as 

Sx,y = j^tr[SFix)Sl,{y)]; s^^y = ^tr[SF{x) S\{y)] . (4.1) 

We will also have to consider quadrupole operators 

1 1 _ _ _ _ 

<lx,y,u,v = —tr[SF{x) Sl,{y) Sf{u) Sl,{v)]; ^x,y,u,v = —tr[SF{x)S%{y)SF{u)S%{v)]. (4.2) 

In the discussion of the total cross section the quadrupole operator appears in 1/A^ suppressed terms, 
however in less inclusive observables |l^ it can contribute already in the leading order. This will be 
the case in our present discussion as well. 

To define the dipole model limit we assume that the initial state of the projectile \P) contains 
only dipoles. We also make two important approximations. First, we consider only the leading N 
part of the evolution. This drastically simplifies the calculations since in this limit there are no 



transitions between different singlet states during scattering. Technically this means |19| that any 
singlet operator [s or q etc.) satisfies a first order differential equation in rapidity. Second, we 
assume the target mean field approximation in which the target averages for dipoles and quadrupoles 
factorize. 

{s{x,y)s{u,v))T = {s{x,y))T {s{u,v))t; (4.3) 
{q{x,y,u,v)s{z,z))T = {q{x,y,u,v))T {s{z, z))t 
{q{x,y,u,v)q{x,y,u,v))T = {q{x,y,u,v))T {q{x,y,u,v))T 

Physically this approximation is tantamount to assuming no correlations between the target fields. 
This is the same approximation as employed in deriving the Kovchegov equation |5|. 

4.1.1 The Hamiltonian H3 

The action of the JIMWLK Hamiltonian on a function F which depends only on s and q in the large 
N limit can be straightforwardly written as: 

jjJlMWLK^g^ q] = [H, + H,] F[s, q] (4.4) 
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with 



oz-s I , 1 " 

as 
'2tt 



{ [Mx,y;z ~l~ ^u,v;z -^x,u,v,y;z] Q: 



x,y,u,v 



-'x,y,u,v\z Sx,v ^y,u 



L 



x,v,u,y;z ^x,y Su,v 



I x,y,u,v,z 

~l~ J-'x,v,u,v;z Qx,y,u,z Sz,v ~l~ Lx^y^x,v;z Qz,y,u,v Sx,z ~l~ Lx^y^u,y;z Qx,z,u,v ^z,y ~l~ -l^u,y,u,v;z Qx,y,z,v Su,z} 

6 



6q{x,y,u,v) 

Here in addition to the Weiszacker-Wilhams kernel 

(x - z)i (y - z)i 



(4.5) 



Kx 



x,y]Z 



[x — z)"^ (y — zy 



we have introduced the dipole kernel 

M — K -i- K — K - 

^"x,y]Z — ^^x,x,z ^^y,y,z ^^x,y,z 

and the "quadrupole kernel" |3^, |ll| 

' {x - z)i {y - z)i 



K 

J-^y.x 



{x - yf 



y,x,z 



L 



x,y,u,v;z 



{x - zf [y - zf 



{x — z)2 {y — z 

{u - z)i _ {v - z)i 
(u — zY {y — zY 



The kernel L can be expressed in terms of the kernel M as 

1 



^y,u;z ~l~ Mx,v;z ■^y,v;z ^x,u;z ] • 



■'x,y,u,v;z — ^ 

For future use we introduce an additional kernel 

{x - z)i {u - z)i 



■^X,Z,Uy 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(x — z)2 [u — z)"^ 

The equation for s is the Kovchegov equation while the equation for q has been derived in [16, p!s| ]. 

The fact that the equation is first order means that the solution for the evolution of any function 
of s and q which satisfies the JIMWLK evolution in the dipole limit is [pl|, w3] 



with sy and qy solutions of the differential equations 



d^s{x,y) 



as 
2^ 



-^'Ix,y\z \^x,z ^z,y 



"^,11} 



(4.11) 



(4.12) 



x,y,z 



and 



a.9 



d^q{x,y,u,v) = — 
I vr 



\J^x,y\z ~l~ M'a'Uiz u,v,y: 



,v,y;zi Hx,y,u,v 



- L 
+ L 



x,y,u,v,z 



x,y,u,v;z ^x,v ^y,u 



Lx,v,u,y;z ^x,y ^u,v ~l~ L'x,v,u,v;z Qx,y,u,z ^z,v 



x,y,x,v;z \lz,y,u,v 



^x,z ~l~ I^x,y,u,y;z Qx,z,u,v ^z,y ~l~ I-'u,y,u,v;z Qx,y,z,v 



(4.13) 
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with the initial conditions 

SY=o{x,y) = s{x,y), qY=o{x,y,u,v) = q{x,y,u,v) . (4.14) 
The quadrupoles and dipoles are not quite independent, in the sense that 

q{x,y,y,v) = s{x,v), q{x,y,u,x) = s{u,y). (4.15) 



It can be checked that eqs. ( 4.12| . 4T^ ) preserve this constraint. That is if the initial condition at 



r] = satisfies eq. ([4.15 ) so does the solution of eqs.( [4.12 , 4.13| ) at any rapidity r]. 



The dipoles s and the quadrupoles q do not exhaust all color singlet operators that enter our 
analysis. We also have to consider quadrupoles of the mixed type 

qZ,v,u = ^ tr[SF{x) 4(y) SHu) sUv)] ■ (4.16) 
These objects appear at the intermediate stages of the calculation even though at the end of com- 



putations (Eq. ( 3.18D ) we always set S = S. The reason this operator appears is that the projectile 



averaged scattering matrix S^'^ in eq. (|3.5|^ ) depends on the product S{x)S^{x). The dipole con- 
tent of the projectile wave function then means that Sg'^ depends on the dipole operator built from 
S{x)S^{x) which is just a special case of q'^'^. 

For the mixed quadrupole q^'^ 



'lx,y,u,v\s=S ~ Qx,y,u,v] Qx,y,y,v\s=S ~ ^x,v 1x,y,y,x\s=S ~ ^ 

ix,x,u,v ~ ^u,v 1 Hx,y,u,u ~ ^x,y ■ 



(4.17) 



The action of the JIMWLK Hamiltonian on a function F which depends on q^'' in addition to s and 
q in the large N limit can be written as: 



HJi^^^^[S]F[s,q,q''] = (i/, + H, + H,s.) F[s,q,q''] (4.18) 



with 



Hqss — —- / {—[Mx^y^z + Kx^y^z] Qx%,u,v + Kx^y-z Sx,y Su,V + 

^ Jx,ii,u,v,z 



' x,y,u,v,z 

+ {Kx^x;z — Kx,y;z) qz,y,u,V^x,z + {Ky^y-z — Kx^y^z) qx,z,u,V^z,y} ^gss(^^ y U v) ' ^^'^^^ 

We also have to derive the dipole limit of the operator Qf [z, 5] S] which enters both the 
evolution Hamiltonian and the expression for the gluon emission vertex Og . Again using the large 
limits to evolve color singlets we obtain 



2 Q'}[Z, S] Q-[z,S] = ^ [ {- {Mx,z,v,2 + My,z,u,,) qZ,u,v + 



I x,y,v,u 

+ {J^y,z,v,z + Mx, z,u,z) ^x,y Su,v ~l" {J^y,z,v,z ~\~ ■^x,z,u,z ^x,z,v,z ■^y,z,u,z) Qz,y,u,z qx,z,z,v 
~\~ {J^x,z,v,z J^y,z,v,z) Sx,z Qz,y,u,v ~^ i-^y,z,u,z ~ ■^x,z,u,z) ^z,y Qx,z,u,v 

+ 

+ {J^x,z,v,z ~ -^XjZjUjz) Sz,v Qx~y,u,z "I" {J^y,z,u,z ~ ■^y,z,v,z)Su,zQx~y,z,v} T~sS ' (4.20) 

'^"x,y,u,v 
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Note a very important property of eq.( [4.20 ): when acting on a function of s, q and g** the 



operator does not produce any additional objects. In principle one might have expected that the 
operator involves higher mixed multipole operators. If that were the case we could not have limited 
ourselves by considering only the set of operators s, q and q^^ but would have to introduce higher 
miltipoles. However as it stands taking into account that T,q^ at initial rapidity depends only on q*** 
we need not include additional operators into consideration at all. 

Our goal is to use the fact that we have to put S = S at the end of the calculation to the utmost 
and thereby to minimize the number of degrees of freedom which we have to deal with. To this end 
we note that the operator Og does not involve derivatives with respect to s, g and s, q but only with 
respect to q^'^. This allows us to simplify the evolution operator by using 

[ {Hs + H-s + H, + Hg) F[s, q- s, q] ]-^^. -^^ = {Hs + H,) F[s, q;s,q] (4.21) 

and substitute s (q) instead of s (q) in expressions for Og and H^sb. 

Since Og depends on derivative with respect to q^^ we cannot set q^^ equal q from the outset. 
We therefore define the "fluctuation" of q'^^ by 

S S 

'ix,y,v,u = 1x,y,v,u + tx,y,v,u J-^ = TT • (4.22) 

As defined above 

tx,x,v,u — ix,y,u,u — . (4.23) 

This follows from Eq. ( [4.17] ). 

At the end of the computation we have to set tx,y,v,u = 0. We now change variables from q"^^ to 
t throughout. Note that Hg as defined in Eq. (4.5) is proportional to 6/6q at fixed q^^. In order to 



correctly account for the shift we have to remember that 

Now we have to express Z in terms of the dipole s, quadrupole q, and t. The Hamiltonian 
becomes 

H3 = Hs + Hg + Ht (4.25) 

with 

Ht — ~ / { [Mx^y^z ~l~ ^'iu,v;z I-'x,u,v,y;z] tx,y,u,v ~l~ I-'x,v,u,v;ztx,y,u,z Sz,v '\' 

J x,y,u,v,z 

~\~ Lix^y^x,v;z tz,y,u,v ^x,z ~\~ I-'x,y,u,y;z tx,z,u,v ^z,y ~l~ Li^ y u^^-z tx^y^z,v Su,z 

Lx,y,u,v;z {tz,y,u,z Sx,v ~l~ tx,z,z,v ^u,y") x ,y ,u,v;zt z ,y ,u,z tx,z,z,v} TT • (^•^^) 

"txjyjUjV 



Thus the four point function t satisfies the classical equation 

d ^ a., 
dr] 

~l~ -^x,y,x,v;z ^z,y,u,v ^x,z ~l~ ■^x,y,u,y;z ^x,z,u,v ^z,y ~l~ -^u,y,u,v;z tx,y,z,v ^u,z 

^x,y,u,v,z {tz,y,u,z Sx,v ~l~ ix,z,z,v ^u,y') Lx^y^u,v;ztz,y,u,z tx,z,z,v} • (4.27) 



'ix,y,u,v — r) / { [■^x,y;z ~l~ Mu,v;z L'x,u,v,y;z] tx,y,u,v ~\~ I-'x,v,u,v;z tx,y,u,z ^z,v ~l" 
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The Hamiltonian Ht describes propagation of t in the background field s. 

The Hamiltonian eq. (|4.26| ) has the following very useful property. Let us introduce the equal 
point limits of the function t: the two-point function (j)x,y = tx,y,y,x, and the three point function 
'4^x,y,v = tx^y,y^v The functions and are not fully independent and their mutual relation is 
precisely that of the dipole s and quadrupole g, ilJx,y,x = (kx,y Also (t)x,x = '^x,x,v = as follows from 
Eq. ( 4.17 ). The evolution of (p can be found by taking the x = f , y = u m. eq.( 4.26| ) and similarly 
for ip. The equations turns out to be decoupled from tx^y^u,v with all four coordinates different 

5r, 4>x,y = ^ / ^x,y;^ + (t>z,y " (t>x,y + (t>z,y (t>x,z] (4.28) 



- Lx,y,y,v;z [i(pz,y + 1) 1px,z,v - 4>z,y Sx,v] | • (4.29) 

The equation for cj) is equivalent to the Kovchegov equation for 1 + (j). We will come back to these 
equations later. 

The coupling of t to s is the origin of the violation of the AGK cutting rule found in [16|. In 



the following we will be expanding the evolution operator expl—H^Y} in powers of t since t has to 
be set to zero at the end of the calculation. This is obviously tantamount to perturbation theory 
in the last "splitting" term in eq. ( 4.26| ). The highest order of the perturbation theory up to which 



we have to expand is determined by the number of emitted gluons. But first we have to rewrite the 
gluon emission operator eq.( [4.20 ) in terms of t. 



4.1.2 The gluon emission operator Og 

The expression for the gluon emission operator is easily obtain, although it is cumbersome to write 
down. First we define for simplicity 

'^■■^ I Jk{z-z) 



Og{k) = -^ / e'^^'--'^ A{z,z) (4.30) 
z vr 



z,z 



The operator A is linear in functional derivative 5/5t. Counting the functional derivative as a power 
the operator contains terms of order minus one, zero and one in t. 

A{z, z) = A-i{z, z) + Aq{z, z) + Ax<yz, z) (4.31) 

with 

A-x = [ J{x,y,u,v;z,z) ^ 



x,y,u,v 5tx^y^u,v 



{ij^y ,z ,v ,z ~^ ■^x,z,u,z') Sx,y Su,v (J^x,z,v,z ~l~ •^y,z,u,z')Qx,y,u,v 

ix,y,u,v 

+ {^fx,Z,V,Z - My^ z,v,z) Sx,z Qz,y,u,v ~l" {J^y,z,u,z J^x,z,u,z) ^z,y Qx,z,u,v 
~l~ {,J^x,z,v,z J^x,z,u,z) Sz^v Qx,y,u,z ~l~ {J^y,z,u,z J^y,z,v,z)Su,zQx,y,z,v 
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-\- {•f^y^z,v,z ~l~ '^X^Z^U,Z J^x^z^v^z ^z ^u^z) Q.Z ,y ^U^Z Q.X ,z ,z^v\ (4.32) 

f s 

Ao = ix,y,u,v D{x,y,v,u;x,y,u,v;z,z) — 

J x,y,v,u;x,y,u,v;z,z 'J'^x,y,u,v 

,u,z) tx,y,u,v 

I x,y,v,u 

~l~ {J^y,z,v,z ^ J^x,z,u,z ■^x,z,v,z ■^y,z,u,z) Qz,y,u,ztx,z,z,v 
~l~ {/^y,z,v,z ~l" J^x,z,u,z ■^x,z,v,z ■^y,z,u,z)Qx,z,z,viz,y,u,z 
~l~ i.J^x,z,v,z ■^y,z,v,2) ^x,z tz,y,u,v ~l~ {•^y,z,u,z ■^x,z,u,z") ^z,y tx,z,u,v 

~l~ i,J^x,z,v,z J^x,z,u,z) Sz^v tx,y,u,z ~l~ {J^y,z,u,z ■^y,z,v,z) Su,z^x,y,z,v\ T7 (4.33) 



^tx,y,u,v 



■^1 = / V[x,y,u,v]z,z) tz^y^u,zix,Z,Z,V 



x,y,u,v 5ix,y,u,v 

{{J^y,z,v,z ~l~ J^x,z,u,z J^x,z,v,z ■^y,z,u,z)tz,y,u,zix,z,z,v\~Fl (4.34) 
x,y,u,v "''X,y,u,v 

The explicit form of D and J is given in Appendix A. It is also convenient for future use to define 
the "vertices" associated with the two point function (l)x,y 

J^{x,y;z,z) = J{x,y,y,x;z,z) (4.35) 

~ {■^y,z,x,z {Sz^z + ^x,y ~ ^x,z ~ ^y,z) ~^ ■^y,z,y,z {Sy^z + ^z,y ~ ^z,z ~ 1) + 

+ ■^x,z,y,z {Sz,z + ^x,y ~ ^z,y ~ ^z,x) + ■^x,z,x,z {^z^x ~^ ^x,z ~ ^z,z ~ ^) } ■ 

We also separate the part of the vertex D which corresponds to the first line of Eq. ( [4.33| ) (the virtual 
term) taken at the points x = v,y = u, x = u, y = v 

D^{x, y; X, y; z, z) = -{Kz,z,x + Kz,z,y) S{x - x) 5{y - y) . (4.36) 

We will also use 

J{x,y,u,v-k) ^ ^ [ e''^'-'^ J{x,y,u,v;z,z) (4.37) 



z,z 



J^{x,y;k) = ^ / e*'^^^-^") Mx,y;z,z) (4.38) 



z,z 

a. 



D{x,y,u,v;x,y,u,v;k) = — - / e^^i^ D{x,y,u,v;x,y,u,v;z,z) (4.39) 

2vr Jz,z 

D^{x,y;x,y;k) = ^ _ e*^^^"^) D^{x,y;x,y; z,z) = - 5{x - x) 5{y - y) . (4.40) 

The last element we need to know in order to proceed is the "matrix element" Oq which charac- 
terizes the initial projectile wave function. 

4.1.3 The matrix element Co 

In the dipole limit, the wave function of the incoming projectile is decomposed in the dipole basis 
and the S - matrix depends only on the dipole operator 

= (4.41) 
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This implies 

^riS^S] = ^r[QZ,yJ = So^^[l + U,y,y,.] = ^^i^ + <t>.,y] (4.42) 

This expression can be expanded in Taylor series in t (Sq'-'^[1] = 1): 



^0^[^ + tx,y,y,x] 



if"*/" (5 1 

^+127;;^ {U d'x,d\t^^,y^,y^,^^- So^^[l + t]It=o. (4.43) 

m=l U=l Ol^^,y„y^,x, J 



For a projectile containing a single dipole at points x, y the expression is very simple 

+ t] = 1 + t^^y^y^^ = 1 + (l)^^y. (4.44) 

Given these ingredients, we can formulate the rules of calculation of multigluon production cross 
sections in terms of diagrams involving the propagation of the field t in the external fields s and q. 

4.2 Diagrammatics and Feynman rules 

In the dipole limit the generating functional Z reads 

Z[j] = 1 + / DsDqDt5{t)W^^[s,q] x 



X e 



/^^ dv[H. + H, + H,-J,o,(k),] ^^ppj^ ^ _ 2^PP[s]) . (4.45) 



The T,qP[s] term contributes only to the total cross section (j = 0) which is not our interest in this 
paper. Since this term does not depend on t it vanishes under the action of Og and hence does not 
contribute to particle production. We therefore omit Sq'-^[s] in what follows. We also omit the unity 
which vanishes upon differentiation with respect to j. 



It is convenient to get rid of Hg and Hg in the evolution operator by rewriting Eq. (4.45) in the 
"interaction picture" with respect to the Hamiltonian Hg + Hg. 

Z[j] = J Dt5{t)e- ^0 dv{Htls^,gr,] - J,0,(k,ls^,g^])j(k,r,)} + t] . (4.46) 

Here s^^ and are "time" dependent external fields, which are solutions of the corresponding equa- 
tions of motion with initial conditions s^=o = s and q-q=o = q- In the spirit of the mean field 
approximation we have omitted the factor Wy^Q which specifies the distribution of these initial 
conditions. This factor can be restored if desired. 

It is now convenient to organize the calculation as perturbation theory in t. The n - gluon 
production cross section is given by the path integral with n insertion of the operator Og{k, [srjj^r;]) 
and thus contains at most n functional derivatives with respect to t. Since t has to be set to zero 
at the end, this means that the rest of the expressions have to be kept only to order t"". Thus we 



treat the homogeneous terms in the Hamiltonian Ht Eq. (4.26) as zeroth order. Their sum generates 
the "free" perturbative evolution, and this evolution has to be solved exactly. Just like with the 
evolution of s and q we see that any functional of t in this framework evolves as 

F^[t] = F^=o[h] (4.47) 
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where tjj satisfies the equation 

Cis f 

df) tx^y^u,v 7^ / { \-^x,y\z ~\~ ^u,v;z Lx^u,v,y;z\ tx,y,u,v ~l~ Lx^v,u,v;z ix,y,u,z ■^z,vi'^ ~^ 

~l~ Lx,y,x,v;z tz,y,u,v ■'^x,z{Y ^) ~l~ -^x,y,u,y;z tx,z,u,v Sz,yi^ ^) ~l~ I-'u,y,u,v;z tx,y,z,v ■'^u,z{^ ^) 
^x,y,u,v;z {tz,y,u,z Sx,viy ^) ~l~ tx,z,z,v Su,y{^ ^))} (4.48) 

Note that since the evolution of s and t starts on the target side while the evolution of t starts at 



the projectile side, the rapidity of the "background field" s in Eq. (4.48) is Y — rj. The same holds 
for the later discussion of the gluon emission operator. This has to be kept in mind although we will 
rarely indicate the rapidity label of s and q explicitly in the rest of this section. 

Eq. ( [4. 481) has to be solved with the initial condition 

tn=o{x,y,u,v) = t{x,y,u,v) . (4.49) 
Formally the solution can be written as 



tr,{x,y,u,v) = / Gr,{x,y,u,v\x,y,u,v) t{x,y,u,v) (4.50) 

J x,y,u,v 



with the propagator 



V 



G,j{x,y,u,v\x,y,u,v) = {x,y,u,v\V exp{— / dri'h[s{Y — ri'),q{Y — ri')]}\x,y,u,v) (4-51) 







where h is the "first quantized" Hamiltonian that can be read off eq.( [4.48 ). In other words the 
propagator G is the solution of the differential equation 

— G{x\x,y,U,v) = J {-[Mx,y-z + Mu,v;z - Lx,u,v,y;z]G{5i\x,y,u,v\) + 

~l~ Lx^v,u,v;z G'(x|x, y, U, z) Sz^v ~\~ I-'x,y,x,v;z G(x|z, y, U, f) Sx^z ~l~ Lx^y^u,y;z G(x[x, Z, U, f) Sz^y ~\~ 

+ Lu,y,u,v;z G{5i\x, y, Z, V) Su,z - Lx,y,u,v;z [Gi5i\z, y, U, Z) Sx,v + ^(xlx, Z, Z, V) Su,y]} (4.52) 

with the initial condition at ?7 = 

Gr,=oi^\x,y,u,v) = dx,xSy,y6u,uSv,v ■ (4.53) 

Here x = {x, y, u, v). 

We will also need the BFKL Pomeron Green's function corresponding to propagation of the two 
point function (j) (now x = (x,y)): 

^G^^^^{5L\x,y) = ^ j Mx,y.,z [G^^^^'^(x|x,.) + G^^^^{5L\z,y) - G^^^'^(x|x, y)] . (4.54) 

The initial conditions for the BFKL propagator: at r/ = 

G™(x|x,y) = (4.55) 
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The last term in eq.( 4.26 ) increases the power of t and thus is to be treated as a perturbation. 
It defines the perturbation vertex 1^2 

6 f a<, ^ . 5 



A(x, X, x) t(x) t(x) 



5t{x) 



x,y,u,v,z 



27r 



Lx,y,u,v\z t 



z,y,u,z tx,z,z,v 



6t 



(4.56) 



x,y,u,v 



We note that when this operator acts on a functional of the two point function 4>x,y only, it 
produces another function of (/) according to the rule (Eq. ( [4.28| )) 



— M, 



x,y,z 



2-K 



which defines a triple (j) vertex 



X^{x,y,z) 



x,y,z 'f'z,y YXyZ 



a 



^x,y 



(4.57) 



(4.58) 



Finally the gluon emission operator defines three types of insertions. The term A^i annihilates 
the propagating field t with amplitude J[s{Y — rj),q{Y — rj)], the term Aq inserts the vertex D[s(Y — 
r]),q{Y — r])] on the propagator of t and Ai splits an existing propagator with the amplitude y [s(y — 
ri),q{Y -T])]. 

This defines the set of Feynman rules for the diagrammatic representation of the "partition 
function" Z eq.( [4.46| ). We denote each factor of t in the expansion Eq. ( 4.43| ) by an empty circle. 



The propagator G of Eq. (4.52) is denoted by a straight line. The propagator carries an arrow 
indicating the direction of "time" Y (Fig. 1). 



G 




Figure 1: The propagator and the vertex A 



The vertex A, Eq. ( 4.56| ) is denoted by a full circle. Every vertex has two propagators entering 
and one propagator leaving (Fig.l) 

Finally we denote insertions of J, D and V of Eq. (4.32),Eq. ( 4.33D ,Eq. (4.34) by a full circle 
with horizontal line. Note that the insertions of J and D are quite complicated, as they depend 
on rapidity dependent fields s{Y — rj) and q{Y — rj). We cumulatively represent the factors s and q 
entering the insertion J by dashed lines (see Appendix A for more details). 

To compute the n gluons inclusive cross section we proceed as follows. 
1. At the projectile rapidity Y draw 1 < m < n open circles corresponding to the expansion of 



Sg'^ Eq. ( [4.43 ) to m-th order (if the projectile is a single dipole there will be only one circle). 



2. Each open circle is a sink for one propagator G Eq. ( fl.52 ). Draw the diagrams of the field 
theory with propagator G and vertex A of Eq. ( [4.56| ) with n insertions of the operators Og{ki) 
an rapidities Y^ + Yi. 
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3. The only diagrams that contribute are the ones where each propagator starts on the source J 
and ends either on the vertex A or on the vertex V or on the open circle. 

4. Finally we sum over all diagrams with 1 < m < n and Ny + Nj + = n, where Ny, Nj and 
Nd are the numbers of insertions of V, J and D respectively. 

In each diagram the number of various vertices are related by 

m + Nx + Nv = Nj (4.59) 
where Nx is the number of vertices A. 

Since there are at most n insertions of the source term J, the perturbation theory in A has to 
be developed at most to order n — 1. 

In the next subsections we explicitly consider the examples of one-, two-, and three-gluon inclusive 
cross section. For simplicity in all these examples we take the projectile as a single dipole (xcyo)- 

4.3 Single inclusive gluon production 

Single gluon production rate is our first example. Since J = 1 and m = 1 there is only one diagram 
(Fig. |2|) which contributes to this observable. Moreover the propagator entering this diagram is 



Yl 



* % 
* \ 
* \ 
* \ 
* % 
* \ 



Figure 2: Single inclusive gluon production 

simply the BFKL propagator. The reason is that the projectile is the functional only of the 
two point function 4'xo,yo- other hand as discussed above, the two point function propagates 

independently of other components of t, Eq. ( [4.28| ). The homogeneous part of Eq. ( 4.28| ) is simply 
the BFKL equation, and so the propagator of (p is the BFKL propagator. 

The analytic expression for the diagram of Fig. 2 is 
da 



dYi dk^ IT 
a. 



z,z,x,y 

ik(z-z) I J^^ 



BFKL 



{xQ,yo\x,y]Y - Yl) J{x,y;z,z;Yi) 



z,y,z 



G 



BFKL 



{xo,yo\x,y;Y - Yi) x 



x,y 



slyiYl 



(4.60) 



This reproduces the result of Kovchegov and Tuchin [21|. 
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4.4 Double inclusive gluon production 



For the double inclusive gluon production there are two diagrams depicted on (Fig. |^ . The diagrams 
are topologically the same as those of |24|, however the propagators in Fig. 3 are not necessarily 
those of the BFKL Pomeron. The diagramm on Fig. 3a is given by 




(a) (b) 
Figure 3: Double inclusive gluon production 



da 



dYi dkf dY2 dk^ 



2 1(a) 



Y-Y^ 



dr] G^^^^^(xo, yo\x, y, v) A^x, y, z) G^^^^(x, z\ui,vi-Y -Yi-r]) 
X G^^^^{z,y\u2,V2;Y -Y2-ri) J^{ui,vi; ki;Yi) J<^(n2, ^2; Ajs; ^2) (4.61) 
where we used Einstein's summation convention for transverse coordinates x,y, etc.. 



The diagram on Fig. ^ is 
da 



dYi dkl dY2 dk 



2m 



G 



BFKL 



(xo, yok, y;Y - Yi) D{x, y, y,x\x, y, u, v; fci; Yi) x 



X G{x,y,u,v\xi,yi,ui,vi;Yi - Y2) J{xi,yi,ui,vi; k2;Y2) . (4.62) 



Note that in this diagram the BFKL propagator of two point function (j)x,y become a propagator 
of the three point function ipx,z,z, '4^z,y,z s.tc. below the insertion of the first gluon emission operator. 
This does not happen in the single gluon inclusive calculation nor in the diagram Fig. |a. The 
propagation of the three point function couples to the "background fields" s and q, and therefore this 



diagram violates the AGK cutting rules as stressed in |lf:]. There is one contribution to this diagram 
however which does not change the nature of (j). This is the insertion of D^{x, y, x, y; k) of eq.( [4.36| ): 

da 



dYi dkf dY2 dk^ 



G 



BFKL 



X G 



BFKL 



{xo,yo\x,y;Y -Yi) D^{x,y;x,y;ki) x 
{x,y\u,v;Yi - Y2) J(n, ^2; ^2) • 



(4.63) 



Nevertheless it does not seem useful to separate this particular contribution as it is not the only 
contribution even in the BFKL limit. 



Our result coincides with that of The detailed correspondence is as follows. The expression 



dT]G 



BFKL 



{xo,yo\x,y;ri) \'>'{x,y,z) G 



BFKL, 



X, z\ui,vi]Y —Yi—ri) G 



BFKL, 



z,y\u2,V2;Y-Y2-r]) 
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is identical with n2 of Ref. Thus the contribution of diagram Fig. 3a reproduces the first term 

in Eq. (32) of that reference. Eq. (|4.63| ) is identical with the very last term in Eq. (32) in [^] , the 
term proportional to M (without the factor S). The expression 



G 



BFKL 



{x,y\u,v\Yi -Y2) J{u,v\k2]Y2) 



is M{x,x,y,Yi;k2,Y2) + (x ^ y) of Ref. |T|1. Finally Eq. ( |462D with Eq. ( ^ ) subtracted is the 
same as the "MS"' terms of Eq. (32) in |16|: the vertex J{ki,Yi) brings in the factor s {S of Ref. 
fl^); the three-point function M of |jl6| equals 



G{x,y,u,v\xi,yi,ui,vi;Yi - Y2) J{xi,yi,ui,vi; k2;Y2) 
whereas the vertex J{k2,Y2) equals d in Eq. (27) of []l6| . 



4.5 Triple gluon inclusive production 

It is straightforward to draw Feynman diagrams which contribute to the triple gluon inclusive pro- 
duction cross section. There are nine distinct diagrams that are depicted in Fig. 4. We do not write 
down analytic expressions for the different contributions: this is a straightforward if cumbersome 
endeavor. We only note that the new element that appears here relative to the single and double 
gluon production is the last diagram which involves insertion of the vertex Ai. Since the vertex Ai 
turns a three point function ip into a generic four point function t, this diagram contains the general 
propagator G below the rapidity Yi. Thus even though the diagrams again are topologically the same 
as in the BFKL Pomeron calculus ||2^, the propagators are those of the generic four point function 
tx,y,u,v rather than BFKL propagators of (px,y 



4.6 Can v^e make sense of it? 

The diagrammatics involved in the calculation of multiple gluon production described in this section 
is simple, however an actual calculation would be quite involved. Single gluon production apart from 
the solution s of the BK equation, requires only the knowledge of the propagator of the two point 
function (px,y which is simply the BFKL propagator and thus is reasonably straightforward. The 
double inclusive gluon requires already the knowledge of the propagator of the three point function 
'4^x,y,v This we do not know in general since ip couples to the "background" field s. For three 
and higher gluon production we would need the knowledge of the propagator of the generic four 
point function t which is even more complicated. In this subsection we motivate a scheme which 
approximates the propagator of i by a known function of the BFKL propagators. 

Let us first examine the evolution of the three-point function ipx,y,v more closely. We rewrite the 
linear part of the equation of motion of ijj here for convenience: 

~\~ Lix^v,y,v;z '4^x,y,z Sz,v ~\~ I-'x,y^x,v\z'4^z,y,v Sx,z Lx^y^y^v\z '4^x,z,v^ • (4.64) 

We are interested in the high energy limit, where the target probed at any rapidity between Yq and 
Y is black. The leading behavior of then is determined by taking the limit s — > 0. We have to be 
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Figure 4: Three- gluon production 



somewhat careful since s{x, z = x) = 1 even in the deep saturation regime and this regulates the UV 
singularity in the virtual term in the equation. Still to get a qualitative idea we can take 5 = but 
use Qs{Y — rj) as the UV cutoff wherever a UV cutoff is required. Taking naively s = in Eq. ( [4.64| ) 
we discover that only the first and the last terms survive. The first, virtual term is negative and UV 
divergent. Assuming that it gives the leading contribution to the evolution we find. 

'4'x,y,v;n ~ exp { - g-r? ln[(x - vf Ql{Y - t?)] } . (4.65) 



This of course is only valid for |x — w| > Q^^iY — rj) since otherwise the last term in Eq. (4.64) gives 
a large contribution which cancels the putative UV divergence. Note that the "source" term - the 
second term in Eq. ( |4.64 ) gives a competing contribution. At large distance |x — w| > Qj^(Y — rj), 



s{x,v) decays very steeply as |33 



s,.,,(y - r?) oc exp { - ^{Y - rf) ln[(x - vf QI{Y - r?)]} (4.66) 

so that close to the projectile rapidity where r] < Y — rj the contribution of this term is subleading. 
Closer to the target rapidity it dominates, still qualitatively the behavior is very similar. We thus 
conclude that while the two-point function (p grows with energy with the BFKL intercept, il^x,y,v at 
large |x — f | > Qj^{Y — rj) decreases with energy very quickly. 
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It is also possible understand the behavior of '4'x,y,v at small values of the separation x — v. For 
X = V, as explained above we have ipx,y,x = 4'x,y When the separation is smaller than Qj^ but 
nonvanishing, we can estimate the behavior of Tp by taking s = 1 in Eq. ( [4.64 ). This is obvious in 
the second term in Eq. ( [4.64D , since it is directly proportional to Sx,v This also gives the leading 
contribution in the third (fourth) term, since the z- integral is peaked at z = v (z = x). Setting 
s = 1 in Eq. ( ^.64| ) we obtain the equation 



a.. 



drj '4^x,y,v — 2 ^ 



1 



L 



x,y,y,v;z 'rz,y 



v,y,v;z 



'4'x,y,z ~\~ I-'x^y^x,v;z '4^z,y,v 



L 



x,y 



,y,v;z '4^x,z,v ^ 



It is easy to check that this equation is solved by 



x,y,v 



^x,y 



'y,v 



(4.67) 



(4.68) 



if 



x,y satisfies the linearized version of Eq. (|4.2S| ) - the BFKL equation. 

In fact the relation Eq. ( [4.68| ) is required by the definition of tx,y,u,v Eq- ( [4.22| ) in the BFKL 
limit, where S{x) = 1 + ia{x) and S{x) = 1 + ia{x) are close to unity. Expanding qx%,y,v to second 
order in a and a and taking the correlator of a to be equal to correlator of a (in accord with the 



discussion leading to Eq. ( 4.21 )) we find the relation Eq. ( 4.68 ) directly from the definition of t. The 
same argument more generally gives 



tx,y,u,v 



x,u ~l~ 



'y,v 



4>y,u) 



(4.69) 



This again can be checked to solve the equation of motion Eq. ( [4.48| ) with s = 1. 

We stress that Eq. ( [4.68 ) and Eq. ( |4.69 ) hold only in the linearized regime for s = 1. Nevertheless 
this means that at small separations, |x — f| < Q^"^; \y — u\ < Q^^ the propagator of the four point 
function t behaves as 



G{xo,yo,uo,vo\x,y,u,v;r]) = 

\ (G^^^^(xo,nojx,u;r?) + G^^^^(yo, t^o|y, ^; r?) 



G 



BFKL 



{xo,VQ\x,v]r]) 



(4.70) 

G^^^'^(yo,no|y,^x;r?)) 



Next let us take a closer look at the gluon emission vertex Eq. (4.32), Eq. ( [4.33| ), Eq. ( 4.34| ). We 
will zoom into one of the terms in the vertex, but the same considerations exactly go through for all 
other terms. Consider one of the terms in Aq 



+ AC, 



Mx 



■^y,z,u,z) Qz,y,u,z ^x,z,z,v 



6tx 



(4.71) 



,y,u,v 



Let us assume for a moment that the coordinates of the field tx,y,u,v entering this vertex are close 
together, so that |x — t;| < (57^; \y — u\ < Q^^ where Qs is defined at the rapidity y — 77 of the 
quadrupole operator q in Eq. ( 4.71 ). Since qz,y,u,z ^ 1 for \z — z\ > Qj^, the coordinates of the 



outgoing field t are also pairwise within the same distance scale ^. The same is true for all the 
other terms in ^q- 
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Thus the vertex Aq itself picks only the region of the phase space in the propagator G, in which 
the propagator can be taken as Eq. ( 4.70| ). 



The vertex Ai is different in this respect, since the coordinates z and z of the outgoing field are 
not necessarily close to each other. This vertex appears only for triple and higher gluon production. 
For these multigluon observables it is therefore necessary to include the large distance suppression 
factor in the propagator itself. As we have discussed above this suppression is the property of 
the linearized evolution equation for t. In lieu of exact solution of this equation one can resort to 
modelling. The following form suggests itself as a reasonable approximation,: 

G{xo,yo, uo, vo\x, y, u, v; rf) = ^Sx,t,(^ - 'n)sy,u{^ - v) ^ (4-72) 
G^^^^(xo, no|x, n; r,) + G^^^^(yo, vo\y, v; r,) - G^^^^(xo, v^\x, v; r,) - G^^^^(yo, noly, u- r?) 

It obviously has correct behavior both at small and large separations and the crossover in the behavior 
occurs on the correct distance scale. 

To summarize, the diagrams for single gluon production contain only the propagator of 4> and 



thus are the same as in the BFKL Pomeron calculus |24]. For double inclusive gluon production 
since the vertex Ai does not appear, the leading contribution is again given by the diagrams where 
the propagators of ip can be taken at s = 1^ 

G^{xQ,yQ,vo\x,y,v]ri) = (4.73) 
\ (G^^^^(xo, yo|x, y; ry) + G^^^^(yo, vo\y, v; rj) - G^^^^(xo, ^;ok, v; rj)) 

In this approximation the propagators again become those of the BFKL Pomeron and the leading 
contribution preserves the AGK cutting rules, modulo the fact that the field s couples directly in the 



vertex D{ki). This conclusion agrees with the argument of [24| 



However for triple gluon production and higher we cannot use BFKL propagators to approximate 
the propagation of ip and t in the external field s. Thus the interaction with s appears not only at 
the gluon emission vertices and the AGK cutting rules are irretrievably lost. 

However, we believe that if the rapidities of all measured gluons are far from the target, a good 
quality approximation for all momenta k except possibly ki ^ Qs{Y — Yi) can be obtained by using 
the Feynman rules with the vertices given in the previous subsections using as the propagator of the 
field t Eq. ( 4.72| ). The calculation then requires only the knowledge of the BFKL propagator, and 
naturally the solution of Kovchegov equation Sx,y 



^There is a slight subtlety here, since the BFKL propagator leads to diffusion in transverse space. Thus one can 
worry about the fact that in Fig. 3b the transverse spread between the points in the propagator entering the vertex 
J{k2) is not necessarily small even if at the vertex D{k\) the spread was within Q7^- If that where the case our 



argument about the reliability of the approximation Eq. (4.73) would be invalid. However the diffusion in the BFKL 



equation is rather slow, so that the growth of radius is slower than the decrease of Qa {Y — rf) as r) grows towards the 
target. It therefore follows that even allowing for diffusion, the transverse spread always stays small enough so that the 
approximation s = 1 can be used for the propagator. 
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5. Discussion 



In this paper we have developed the formahsm to calculate multigluon inclusive cross section at high 
energy for arbitrary number of gluons. The formulae in the full unabridged JIMWLK formulation 
are complicated and at the moment it is not clear how to develop a feasible calculational approach 
in this regime. We have also discussed the dipole model limit of the JIMWLK approach. In this 
case the expressions simplify dramatically. We have been able to develop a simple diagrammatic 
technique based on Feynman rules for calculations in this model. We have shown that for the double 
gluon inclusive production one recovers the results of |16|. We note that the formalism can be 



extended easily to other observables. In Appendix B for example we consider gluon emission with 
fixed momentum transfer of the valence part of the projectile wave function. 

We confirm that the AGK cutting rules are violated for multigluon production. The violation 
is due to the appearance of non Pomeron exchanges even within the dipole model. The calcula- 
tion involves propagators of three- and four point "multipoles" propagating in the external field s. 
However we presented an argument that for the double inclusive gluon production the leading high 
energy behavior is indeed given by the diagrams preserving the AGK cutting rules in accordance 



with conclusion of |24]. On the other hand for three and higher gluons the AGK rules are genuinely 
violated. 

Perhaps the most remarkable result of this paper is that only finite number of multipoles (namely 
multipoles which depend on two, three and four transverse coordinates) is needed for the calculation 
of any number of emitted gluons. 

The propagators of these multipoles except for the dipole are not known analytically, but the 
problem is still on an entirely different level of complexity than the full JIMWLK. We have suggested 
an ansatz for the propagators based on its known UV and IR behavior which we believe should be a 
very good approximation for calculating the cross section. 

It is interesting to note that although the vertex in the Hamitonian Eq. ( [4.26| ) contributes only 
perturbatively to a fixed order when the number of gluons is fixed, one could in principle calculate 
other observables where the full evolution of t will enter the fray. An example of such an observable 
could be average multiplicity. Since the evolution of t is nonlinear and undoubtedly leads to saturation 
of growth of t this suggests that the average multiplicity must be growing much slower with energy 
than in the pure BFKL calculation. This slowdown is of course affected by the saturation of s but 
mostly by the saturation of t. It would be interesting to study this question in detail. In particular it 
is possible that classical solutions of the evolution equation for t could be directly related to physical 
observables, much like classical solutions of field equations are related to particle production in QCD 
l2c 
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A. Appendix: Vertices 

In this Appendix we collect explicit expressions for the vertices J and D which enter the gluon 
emission vertex eqs.( 4.32| , 4.33 ). 



Ji'^iUj'^j'^ • -2) -2) — {-^yjZjVjZ ~^ ■^x,z,u,z) Sx,y Su,v {/^x,z,v,z ~\~ ■^y,z,u,z)Qx,y,u,v 

~l~ ij^x,z,v,z ■^y,z,v,z")Sx,zQz,y,u,v ~l~ {J^y,z,u,z J^x,z,u,z) Sz^y Qx,z,u,v 
~1~ ij^x,z,v,z ■^x,z,u,z) Sz,v Qx,y,u,z ~l~ {J^y,z,u,z J^y,z,v,z) ^u,z Qx,y,z,v 
~l~ ij^y,z,v,z ~\~ ^x,z,u^z J^x^z^v^z J^y,z,u,2) Qz,y,u,z Qx,z,z,v 

(A.l) 

For the two point function (p: 

J(j>{x,y', Z, z) = My^z,x,z {Sz^z ~^ ^x,y ~ ^x,z ~ ^y,z) ~^ ■^x,z,y,z {s + ^x,y ~ ^z,y ~ ^z,x) (^-2) 

2 I 2 
y,z ^ ■^z,y 



+ ■^y,z,y,z {^y,z + ^z,y ^z,z ^) + ■^x,z,x,z {^z x + ^x,z ^z,z ^) 



* * 

♦ ♦ 
♦ 




♦♦ /\ /» ■ '» 

♦ 



q q q q 



Figure 5: The "source" term J in the gluon emission vertex. 
The vertex D is given by 

D{x, y, U, V, X, y, U, v; Z, z) = {Mx,z,v,z + ■f^y,z,u,z)Sx,xSy,yK,u^v,v (A.3) 
~l~ {J^x,z,v,z J^y,z,v,z)Sx,z^x,z^y,y^u,u^v,v ~l~ ij^y,z,u,z J^x,z,u,z") Sz,y ^x,x ^y,z ^u,u ^v,v 
~l~ {J^x,z,v,z J^x,z,u,z) Sz^v ^x,x ^y,y ^u,u ^v,z ~^J^y,z,u,z J^y,z,v,z) Su,z ^x,x ^y,y ^u,z ^v,v 
~l~ (J^y,z,v,z ~l~ ■^x,z,u,z J^x,z,v,z J^y,z,u,z)iSlz,y,u,z5x,x^y,z^u,z^v,v ~l~ (lx,z,z,v ^x,z Sy,y ^u,y ^v,z) 

Finally the vertex V is given by 

V{x,y,u,v;z,z) = {J\fy^z,v,z + J^x,z,u,z - ■N'x,z,v,z - J^y,z,u,z) (A.4) 



The vertex Eq. ( [A.2| ) coincides with s in eq. (20) of [16|. 



B. Appendix: Multi-gluon production with momentum transfer 

It is straightforward to extend the formalism discussed in this paper to the case when in addition to 
emitted gluons, the transverse momentum transfer of the valence part of the projectile wave function 

''Note that the notations in this paper differ from those in jl6| so that the letter s denotes a different quantity than 
in @ 
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is fixed. The only modification is ttiat of tlie matrix element Oo[S', S]. To impose momentum transfer 
p we take [||] 

Oo[S,S]p = j SB<fae'^''^^P[S\x)S{x + a)] (B.l) 
Here B stands for impact parameter of the projectile. In the dipole limit this becomes 

Oo[S,S]p = j SBcfae''P^i:PP[l + t^,y,y+a,x+a] (B.2) 

One important difference from the formulae discussed above is that the four-point function t enters 
already on the level of the single inclusive gluon cross section. 

For a single dipole projectile {B = (xq + yo)/2) the diagrams are topologically identical to those 
given in Section 4. The only difference is that the "sink" on which the propagators end is given by 
Eq. ( |B.2| ). The expression for the single inclusive gluon production is 

dN 



dYi dk"^ d^p 



j d'^xod^ae'^"- G{xo,yo,yo + a,xo + a\x,y,u,v;Y - Yi) J{x,y,u,v; k,Yi) (B.3) 



We can also consider in this framework the case when the projectile consists of a single quark. 
The observable then corresponds to the process q + A ^ q{p) + X + g{ki; Yi) + g'(fc2; I2) + •••• To 
this process we take as the projectile single dipole with the transverse coordinate of the antiquark at 
infinity: 

Oo[S,S]p = J d^xod^ae'P'^il + (B.4) 
The result for the single inclusive gluon production is 



dN 



dYi dk'^ d^p 



j d^xod'^ae'^" G{xo,yo,yo + a,xo + a\x,y,u,v) J{x,y,u,v;k,Yi)\yg^oo (B.5) 



The propagator G{xQ,yQ, y^ + a, + a|x, y, n, v) |yo-^oo diverges in the infrared and has to be regular- 
ized. This is of course the reflection of the fact that the cross section for a color nonsinglet projectile 
is logarithmically IR divergent. If the gluon is emitted at the same rapidity as the valence quark the 
expression becomes 



da 



dYi dk"^ d^p 



j d^Bd^ae'P^ J{xo,yo,yo + a,xo + a;k,Y)\y^^oo (B.6) 



which agrees with the result of |17|. 
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